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Abstract 



The interactions responsible for Feshbach resonances in binary atom col- 



^ , lisions produce a condensate of molecules in the atomic BEC-system. We 

b : 

discuss the 'ground state' of the condensate system and illustrate that its 
properties differ quantitatively and qualitatively from the properties of a sin- 
gle condensate with atom-atom interactions described by the effective scat- 
tering length that characterizes the binary atomic collisions. We show, for 
example, that the tunneling of atom pairs between the atomic and molecular 
condensates lowers the energy and can bind the condensate system, giving it 

the liquid-like property of a self-determined density. 
PACS numbers(s):03.75.Fi, 05.30.Jp, 32.80Pj, 67.90,+z 
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The recently observed atomic Bose-Einstein condensates |IJ are dilute superfluid quan- 
tum gas systems and their many-body properties are highly interesting. In contrast to the 
situation in the traditional (Helium) superfluids, the strength of the inter-particle interac- 
tions in the atomic condensates can vary over a wide range of values. In particular, the 
scattering length that characterizes the atomic interactions can be negative, corresponding 
to an effective inter-atomic attraction. A condensate of such bosons collapses in the ab- 
sence of an external potential, but a trapping potential may hold a small negative scattering 
length condensate in a metastable state || that can tunnel ||. Moreover, the remark that 
the scattering length may be varied continuously for binary collisions seemed to promise a 
'tunable' interaction strength. The proposed schemes Q to change the scattering length 
involve a resonant scattering process with an intermediate bound molecular state of the 
interacting atoms. In the scheme that involves a Feshbach resonance ||, recently achieved 
experimentally in an atomic BEC-system |J, the molecular intermediate state energy is 
tuned near-resonance (i.e. near the continuum of the incident atom channel) by varying an 
external magnetic field. 

In this letter, we point out that the interactions responsible for the Feshbach resonance in 
the binary atom collisions produce an additional condensate of molecules. Near-resonance, 
the properties of the hybrid atomic/molecular condensate differs qualitatively from the prop- 
erties of a single condensate. As a particularly striking illustration, we will show that the 
tunneling of atomic pairs between the atomic and molecular condensates can produce a 
bound many-body state with the liquid-like property of a self-determined density. 

The binary atom Feshbach resonances studied by varying a strong external magnetic 
field in an alkali-atom trap are hyperfine-induced spin-flip processes that bring the colliding 
atoms to a bound molecular state of different electron spin. The interaction of two alkali 
atoms, their spins aligned by the magnetic field, is, to lowest order, described by the triplet 
molecular potential. If the magnetic field is tuned near-resonance, the hyperfine interaction 
of the electronic (s) and nuclear (i) spins, V^f ~ s • i, can simultaneously flip the electronic 
and nuclear spins of one of the colliding atoms, bringing the molecular system to a bound 
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vibrational eigenstate of a singlet potential of nearly the same energy as the continuum of 
the incident atom channel. For the ultra-cold condensate systems, the matrix element for 
the transition of the initial state of atoms 1 and 2 in spin state \Si n ) to the bound molecular 
state of spin \S m ), can be taken to be independent of the momenta of the atoms. This 
matrix element, a, is equal to a = (S- m \V h f(l) + V h f(2)\S m ) x / d 3 x ^(x)^(x), where x 
denotes the relative position of atoms 1 and 2 and ip m and ip represent the bound state and 
continuum wave functions of the initial and intermediate molecular states. The continuum 
wave function is normalized by requiring ip to go over into the partial s-wave component of 
a plane wave at large separation. The hyperfine-induced spin-flips are then described by the 
following Feshbach resonance contribution to the Hamiltonian: 

H FR = aJ d 3 r ft m (r)ft(r)ft(r) + h.c, (1) 

where (ft) and (ft) represent the usual annihilation and creation field operators of the atoms 
(ip a ,ft a ) and the molecules in the bound vibrational state (i/j m ,ft m ). 

For the binary atom system, H FR gives a resonant contribution to the atom-atom inter- 
action strength A a (where A a is proportional to the scattering length a s , \ a = A7ih 2 a s /m): 
A e ff = A a — 2a 2 / e, where e is the energy of the intermediate molecular state relative to the 
continuum of the incident atoms. The value of e, which we shall refer to as the 'detuning', 
linearly depends on the magnetic field H. To see that, we note that the energy difference 
A of the singlet and triplet continuum levels, which is the energy gain experienced by a 
single atom in the spin flip process, is proportional to H: A = fi'H. The detuning is the 
difference of the A and the the binding energy E b of the resonant singlet molecule, so that 
e = ji'(H — H ), where H = E b //j,' is the resonant magnetic field. The expression for A c ff is 
most easily obtained by determining the energy shift of two atoms with the same momen- 
tum, confined to a volume Q. In first-order perturbation, the atom-atom interaction term, 
[X a /2]ftftftft, contributes A a /f2 to the energy shift AE . In second-order perturbation 
theory, Hp R couples the initial two-atom state |ini) of energy E in[ to the intermediate molec- 
ular state |int) of energy E int , contributing \(mi\H FR \'mt)\ 2 /[E- mi — E- mt ] = [— 2a 2 /e]/VL . The 

3 



overall shift is AE = \ c ff/Q, where A e ff = A a — 2a 2 /e 0. 

In the binary collision, the intermediate state is virtual since it does not satisfy energy 
conservation. In the many-body system, two particle energy conservation considerations do 
not apply and we argue that, under the appropriate conditions, Hfr produces a condensate 
of molecules. If all atoms are in the same center-of-mass state, then the molecules created 
in the Feshbach resonance also occupy a single center-of-mass state, the characteristic of 
a condensate |§. We can describe the dynamics of the atomic and molecular condensate 
fields, <f) a and m , with an effective Hamiltonian-density ||T0| : 
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where the r-dependence of the 0-fields is understood, and where A a , A m and A denote the 
strengths of the atom-atom, molecule-molecule, and atom-molecule interactions. Having 
produced the molecular condensate, Hfr describes tunneling of atom pairs between the 
atomic and molecular condensates. In the tunneling, the total number of atomic particles is 
conserved if we count each molecule as two atomic particles. We find it useful to introduce 



the atomic particle density n(r) = n a (r) + 2n m (r), where n a 



I 2 and n r 



are 



the atomic and molecular densities. Note that the effective Hamiltonian depends on the 
phase difference of 2 and m . In describing the static situation, choosing a to be real 
and positive, and noting that the if^R-contribution is minimal when m is real, it follows 
from Eq.(Q) that we can describe the statics of a homogeneous condensate system with the 
following energy density: 



u 



A a 2 i ^rn 2 , \ , n 

—n a + — n m + Xn m n a + en m + 2an a (, 



To minimize the energy, we introduce a variable x, which is a scaled molecular field, 
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xyn/2. The square of x then represents the fraction of atoms that have been converted to 
molecules, n a = n(l — x 2 ) and the energy per particle, e = u/n, reduces to 

e = ^n + aV2nx + {[-A + \/2}n + e/2}x 2 - aV2^x 3 + [A /2 - A m /4 - A/2]nx 4 . (4) 



Minimizing e with respect to the molecular field parameter x within the interval (— 1,+1) 
determines the equilibrium value of the molecular field for a given atomic particle density n. 
Taking a to be positive [TTJ , we reduce the problem to the minimization of the dimensionless 
function f(x) = (u — /ay/2n, 

f(x) = x + e'x 2 - x 3 + fix 4 . (5) 

Here e' = {e + n[-2X a + A]} /2« v / 2n, and (3 = n[X a + A m /4 - X]/2aV2n. 

A limit of particular interest is that of off-resonant detuning: e >> 2aV2n, and e 
exceeding the single particle interaction energies, e >> X a n, X m n, Xn. Then f(x) ~i + e'x 2 , 
with a minimum at x = — (l/2e'), giving a molecular density n m pa n 2 [a/e] 2 . The ground 
state energy per particle is equal to u pa nX e s/2, where A e ff is the effective interaction 
strength in the binary collison picture, A e // = A a — 2a 2 /e. Thus the off-resonant limit 
(e >> 0) gives the same result as the effective scattering length description. 

In contrast, as e is lowered near and 'below' resonance (i.e. e —>■ or e < 0), the 
condensate system does not exhibit resonance features. Instead, as the detuning is lowered 
in a condensate of fixed density 0, the fraction of molecules increases until at e = — 2a^/2n + 
n[X — X m /2] jnfl , all atoms are converted to molecules. 

The Feshbach-resonance parameter a, proportional to a molecular overlap matrix ele- 
ment, can take a wide range of values. Calculations for the Feshbach resonances reported 
at MIT, indicate that a^/2n was of the same order of magnitude as X a n. To be definite, 
we will dicuss a system of density n ~ 10 14 — 10 15 cm~ 3 with a\/2n = 5X a n Q , roughly 
corresponding to the Feshbach resonance of largest width detected in the MIT experiments 
0. If a\^2n significantly exceeds the interaction energies X a n, X m n, Xn, then << 1 and 
the fourth order term /3x A of the /-function of Eq.(|^) may be neglected (except when the 
system approaches the purely molecular condensate). Minimizing f(x) pa x + e'x 2 — x 3 , we 
obtain a simple analytical expression for the molecular field paramater x: 

x=[e' - Ve l2 ~T3]/3 . (6) 
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With this result, we obtain analytical expressions for such quantities as the chemical po- 
tential, fi = (du/dn), pressure, P = n 2 de/dn, sound velocity etc... For example, in this 
approximation, (3 = n[X a + A m /4 — \]/2a^/2n w 0, we obtain 



In Fig. (1) we show this pressure for fixed density n = n , and A = 2A a and A m = 4A a 
(implying /3 = 0, so that Eq.(§) is exact), as a function of the detuning. 

The negative pressure near zero detuning indicates that the condensate system increases 
its density. We consider a condensate that is confined by fixed 'walls' to a well-defined volume 
that is large enough to describe the condensate by the above theory (which is an infinite 
system treatment). Lowering the detuning to a value of negative pressure, causes the system 
to occupy a volume less than that allowed by the walls. In the effective scattering length 
description, which gives negative pressure when e < 2a 2 /X a , that volume is zero: the system 
collapses. In contrast, the hybrid atomic/molecular condensate takes on a finite volume, 
determined by the density that minimizes the energy e. In Fig. (2) we plot e as a function 
of the relative atomic particle density n s = n/n for different values of the scaled detuning 
e s = e/ay/2/no for the same parameters as Fig.(l). At e s = 5, for which the effective scattering 
length description predicts zero pressure, the condensate system still gives a positive pressure 
and e is a monotonically increasing function of the density. Without the confining walls the 
condensate expands indefinitely as a gas. At e s = 4.5, the pressure is still positive at n = n Q , 
but the energy is negative at small densities, implying that the many-body system is bound. 
The binding stems from the atom pair tunneling contribution 2an a (p m to the energy density. 
At e s = 4.5, the energy e reaches a minimum at density less than n . Thus, without walls, 
the condensate expands until its density reaches the value that minimizes e. Since that 
density requires a volume greater than that allowed by the walls, a confined condensate 
keeps its density. At e = 3.8, on the other hand, the pressure is negative at n = no and the 
density that minimizes the energy exceeds no: the condensate spontaneously decreases its 
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volume to a fraction of the volume of confinement. The condensate finds its own volume 
and displays the liquid-like property of a self-determined density. Interestingly the value of 
this density (or, alternatively, the volume) is determined by the detuning and can be varied 
by changing the magnetic field. 

In Fig. (3), we show the trajectory that the condensate follows in the plane of the detuning 
and density, as the detuning is lowered adiabatically through the Feshbach resonance. The 
condensate is gaseous and occupies all of the allowed volume - its density remains equal to 
no until its pressure vanishes. Lowering the detuning further, the condensate contracts and 
subsequently expands along the curve of zero pressure. Note that for a range of negative 
detuning values, we find two densities at which the pressure vanishes. The lowest density 
at which this happens actually corresponds to a maximim in the density, rather than a 
minimum, and the system could not remain on this part of the curve. Before the P = 
curve 'turns back', the density n discontinuously changes to its initial value, corresponding 
to an abrupt expansion of the condensate to fill the maximum allowed volume. For an 
atomic-trap with a smooth trapping potential, instead of 'walls', we expect that a similar 
contraction and subsequent expansion. 

The extent to which these very interesting phenomena can be observed crucially depends 
on the relevant time scales. Since the Feshbach resonances create molecules in high vibra- 
tional states, the molecular lifetime, r m , is significantly shorter than the atomic lifetime. The 
decay rate, r" 1 , is proportional to the atomic and molecular densities, r m : = c ma n a + c mm n m , 
where the c-coefficients are the rate coefficients due to molecule-atom (c ma ) and molecule- 
molecule (c mm ) collisions. Estimates based on calculations for hydrogen give molecular 
rate coefficients ~ 10 9 — 10 11 cm 3 /sec | 13| . Using the lower estimate, a condensate of pure 



molecules of density 10 14 cm~ 3 lives for 10 _3 sec. With the relevant time scale for the con- 
densate dynamics h/Xn ~ 10 _5 seconds, it might be possible to adiabatically tune e near 
resonance, although these estimates clearly underline the need for calculations of these 
quantities. Off-resonance (e > 0), only a small fraction of atoms is converted to molecules, 
and, although the molecules decay by colliding with the remaining atoms, the condensate 



is constantly being replenished by tunneling from the atomic condensate. In this limit, 
n m ~ n 2 a \a/e\ 2 and particle loss may be described by a rate coefficient that is the usual 
atomic rate coefficient, c aa , enhanced with the rate at which atom pairs which have been 
converted to molecules are being lost to molecule-atom collisions: dn/dt ~ — c e //n 2 , with 
c e ff = c aa + 2c am n\a / 1\ 2 . 

In conclusion, we have pointed out that the Feshbach resonance interactions can produce 
a condensate of molecules in an atomic BEC system. Off- resonance, many properties of 
the condensate system are described correctly in the effective scattering length description. 
However, this description does not predict the small molecular condensate, which could 
be kept for times exceeding the lifetime of a pure molecular condensate. Near-resonance, 
the tunneling between the atomic and molecular condensates can dramatically alter the 
properties of the condensate systems, a fact that we have illustrated by discussing the 
appearance of a condensate phase with self-determined density. 
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Figure Captions 

Fig.l : Scaled pressure P s = P/[\ a nl/2] as a function of the scaled detuning e s = e/[a^2n \. 
The dashed curve shows the effective scattering length prediction. 

Fig.2 : Plot of the energy per particle, e s = e/[\ a n /2], as a function of the scaled density 
n s = n/rio for the detunings shown on the figure. Negative e s implies a bound many-body 
state. 

Fig. 3 : Plot illustrating the density variation of a condensate, confined by fixed walls, as the 
detuning e s = e/[a v / 2no], is adiabatically lowered through the Feshbach resonance. For the 
lowest initial density, the system contracts and expands along the curve of zero pressure, 
shown in dashed line. Before returning to its initial density, the system undergoes an abrupt 
change in the density. Similarly, the system of medium initial density contracts and expands 
but does not change its density abruptly. The density of the system of highest initial density 
remains constant. 
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